Enhancement of Tc by disorder in underdoped iron pnictide superconductors 
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We analyze how disorder aifects the transition temperature Tc of the s^ superconducting state in 
the iron pnictides. The conventional wisdom is that T^ should rapidly decrease with increasing inter- 
band non-magnetic impurity scattering, but we show that this behavior holds only in the overdoped 
region of the phase diagram. In the underdoped regime, where superconductivity emerges from a 
pre-existing magnetic state, disorder gives rise to two competing effects: breaking of the Cooper 
pairs, which tends to reduce Tc, and suppression of the itinerant magnetic order, which tends to 
bring Tc up. We show that for a wide range of parameters the second effect wins, i.e. in the 
coexistence state Tc can increase with disorder. Our results provide an explanation for several 
recent experimental findings and lend additional support to s^~-pairing in the iron pnictides. 
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Introduction. The symmetry of the superconducting 
state of the iron-based superconductors (FeSCs) is still a 
subject of intense debate^. Photoemission experiments 
on moderately-doped FeSCs show quite convincinglyi^ 
that the pairing state is s— wave, i.e. fully gapped. 
However, since the FeSCs are multi-band systems, the 
s— wave superconducting state (SC) can have either s'^'^ 
symmetry, if the gaps on different Fermi-surface pock- 
ets have the same sign, or s"* symmetry, if the gaps on 
different pockets have opposite signs^. The s"' state 
emerges due to a repulsive inter-band interaction en- 
hanced by spin-fluctuations - this is the key element 
in the theories of a magnetic pairing mechanism in the 
FeSCs^. On the other hand, the s++ state emerges if the 
inter-band interaction is attractive and is enhanced by 
orbital fluctuations!^ 

A seemingly straightforward way to distinguish be- 
tween s"' and s^"*" pairing symmetries is their responses 
to impurity scattering. Both SC states are nearly unaf- 
fected by intra-band scattering. However, while inter- 
band scattering is harmless to the 5"*"+ state, it is pair- 
breaking to the s^ state, leading to a suppression of 
Tci^ To verify this experimentally, one has to choose a 
dopant that acts predominantly as a non-magnetic impu- 
rity scatterer. In the FeSCs, this is not a trivial task, as 
many transition-metal dopants significantly change the 
carrier concentratiouii 

One direction explored by many groups was to substi- 
tute Zn for Fe. Early data on LaFeAs(0i_2:F2:) showed 
that Tc weakly depends on the Zn concentration^, and 
were interpreted as an evidence in favor of an 5"*"+ state. 
Subsequent studies^iiS, however, found that the effect of 
Zn substitution depends on the doping level x: while in 
the overdoped regime Tc displays a sharp decrease, in 
agreement with what is expected for an s^ state, at 
optimal doping Tc remains virtually the same. More sur- 
prisingly, in the underdoped regime of LaFeAs(0i_a;F2,), 
Tc increases with Zn concentration^, an observation that 
is puzzling not only for an s'' SC state, but even 
for a conventional 5+"*" state. A similar increase of 
Tc with disorder was found in the underdoped mate- 



rial Ba(Fei_^Coa;)2As2 with substitution of Cu for Feii. 
That the dopant Cu atoms act as impurity-scatterers fol- 
lows from both band structure calculations^ and neutron 
scattering experiments^^ Intriguingly, measurements in 
the same Ba(Fei_^Co2,)2As2 materials, but with Zn re- 
placing Fe, showed that Tc decreases with increasing Zn 
concentration even in the underdoped regioniS, although 
Tc decreases faster in the overdoped region. 

In this paper, we describe the effect of disorder on the 
superconducting transition temperature Tc of the s'' 
state in underdoped samples, when superconductivity de- 
velops in the presence of SDW order. We argue that 
the conventional wisdom that Tc decreases with increas- 
ing impurity concentration does not work in the under- 
doped region. Indeed, even in a clean system, the rea- 
son why Tc goes down deep in the underdoped region 
is because SDW order competes with superconductivity. 
As doping increases, the SDW order becomes weaker, 
and Tc increases J^ii^ When disorder is added at a fixed 
doping concentration, it influences SDW and SC orders 
differently: while both intra-band and inter-band im- 
purity scattering weaken SDW,— only inter-band scat- 
tering is pair-breaking for s^ superconductivityi^iii As 
a result, disorder without inter-band scattering compo- 
nent does not directly affect SC pairing, but weakens 
the SDW order, leading to an increase in Tc- The sit- 
uation is more complicated when both intra-band and 
inter-band scattering components are present. In this 
situation, impurity scattering affects Tc both directly, 
via pair-breaking, and indirectly, via the suppression of 
SDW. Therefore, the two effects push Tc in opposite di- 
rections, and whether Tc increases or decreases with in- 
creasing impurity concentration depends on the interplay 
between the system parameters. 

Our key results are summarized in Fig. [l] where we 
compare the phase diagrams of FeSCs with and with- 
out impurities for two representative values of the SDW 
and SC couplings and for on-site impurity potential (i.e. 
equal intra-band and inter-band impurity scatterings). 
For one set of parameters. Fig. [TJa). Tc increases with 
increasing disorder in the underdoped region. This be- 
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Figure 1: (Color online) Temperatures of normal (N)-to- 
SDW, SDW-to-SC and N-to-SC transitions as functions of 
doping for the clean and dirty cases (dashed and solid lines, 
respectively) for two sets of system parameters. In the un- 
derdoped region, where SC emerges from a pre-existing SDW 
phase, Tc for s^ pairing increases with the concentration 
of non-magnetic impurities for one set of parameters (panel 
a) and weakly decreases for the other set (panel b). These 
two behaviors are consistent with the data in Refs. [gUllll 
and Ref. [13 |. respectively. Temperatures and So are mea- 
sured in units of Tc.o, which is the SC transition tempera- 
ture at perfect nesting and without SDW (for pure SDW, 
the corresponding temperature is Tn,o)- We used in (a): 
Tm,o/Tc,o = 2, S2/ {2ttTc.o) = 0.4, and impurity-scattering 
amplitudes To = T^ = 0.006 {2tvTc,o); and in (b): Tn,o/Tc,o = 
4, S2/ {2-kTc.o) = 0.8, and To = T^ = 0.012 (27rrc,o). ' 



havior provides an explanation for the experimental re- 
sults of Refs. J9|jll|. For the other set of parameters. 
Fig. [TJb) , Tc decreases with disorder in the underdoped 
region, but with a smaller rate than in the overdoped re- 
gion. This behavior is consistent with the experiments in 
Ref. [11. 

We also analyze the dependence of Tc on the ratio of 
intra-band and inter-band scattering rates (see Fig. [2]). 
For some parameters, as those in Fig. [T]Ja), Tc increases 
for any ratio of intra-band and inter-band impurity- 
scattering. For other parameters, as those in Fig. [Ijb), 
Tc increases when intra-band scattering dominates and 
decreases when inter-band scattering dominates. Even 
in the latter case, the rate at which Tc decreases with 
the strength of disorder is smaller than in a pure s^ 
superconductor . 

The model We follow earlier work a^^'-'^^i^^ and con- 
sider a minimal two-band model for the interplay between 
itinerant SDW and s^ SC. In particular, we consider a 
circular hole pocket at the center of the Fe-only Brillouin 
zone, and an elliptical electron pocket displaced from the 
center by Q = (7r,0) (or (0,7r)). The non-interacting 
fermionic Hamiltonian is 
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Figure 2: (Color online) ATc = T^""'^ - Tf °='" (in units 
of T^^'^^^) as function of the ratio between intra-band and 
inter-band impurity scattering amplitudes (Fo and r,r, re- 
spectively). Line SDWa is for the parameters of Fig. 1(a) 
at a fixed 5o / {2-nTc,o) — 0.2 and line SDWb is for the pa- 
rameters of Fig. 1(b), at a fixed 5o/(27rrc,o) = 0.6. For 
both curves the system is in the coexistence region and 
(Fo + F^)/(27rT=''^='") = 0.01. The lower curve is for a pure 
superconductor, M = 0. In the coexistence region, ATc is 
definitely positive when F.^ is small. For Fo ~ F.^, the behav- 
ior of ATc is a result of the competition between the direct 
pair-breaking effect of impurities, which tends to reduce Tc, 
and the suppression of the SDW order parameter, which tends 
to increase Tc- 



where the operators c (/) refer to electrons near 
the hole (electron) pocket and the band disper- 
sions are given by ei.k = — fc^/ {2m) + fih and 
e2,k = kl/ {2mx) + ky/ {2my) — fie- For small el- 
lipticity \mx — m.yl <^ m, the latter can be conve- 
niently parametrized by S,k = (£i,k ~ £2,k) /2 ~ v^^ • 
(k-kF) and 5k = (ei,k + e2,k) /2 = (yU/i - /ie) /2 + 
(A:^ — kT) {nix — "niy) /4m « ^o + <^2 cos (p. We consider 
the interactions between the low-energy fermions in the 
SDW (particle-hole) and SC (particle-particle) channels, 
as well as their interaction with non-magnetic impu- 
rities. We first introduce the SDW order parameter 

M oc ^[^ ^cr ( Cko-Zkcr ) and reduce the four-fermion SDW 
interaction to 
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where M is obtained self-consistently (see Eq. ([S]) 

below). Introducing the Nambu operators ^^jj^. = 

I C]|^..|,, c_ki, /kt5 /-kj. ) ■ the bare Green's function is ex- 
pressed as 
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where fi are Pauli matrices in Nambu space. In the 
Born approximation, impurity scattering gives rise to the 
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Figure 3: (Color online) (a) SDW transition temperature (in 
units of Tjv.o) as function of doping So (in units of 2ttTn,o) for 
the clean (7 = 0) and dirty cases {■j/2nTM.o = 0.06 and 0.1). 
We set 82 = 0, but the behavior at a finite 82 is quite sim- 
ilar. The dashed lines denote metastable solutions of the 
SDW gap equations, characteristic of the first-order char- 
acter of the transition. In the dirty case, the second-order 
SDW transition extends to T = above a certain thresh- 
old of scattering amplitude, (b) Magnetization M (in units 
of 2-kTm,o) as function of temperature for 5o/ {2ttTn,o) = 0.1 
and 82/ {2-rTm,q) = 0.2 (same as in Fig. [Ha)). The dashed 
and solid lines are for the clean and dirty cases, respectively 
(in the dirty case 'y/{2nTN,o) = 0.006). 



self-energy correction E 
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is the density of impurities, Q 



k^kC/k-q, where 
= Qq^ — S is the 



renormahzed Green's function, and U is the impurity po- 
tential, which we decompose into an intra-band contribu- 
tion uq and an inter-band contribution m^: 



t/ = 



UqTz Ut,Tz 
U-kTz UqTz 



(4) 



To find Q, we write it in the same form as Eq. ([3]), 
but with renormalized parameters d;„, Sq, and complex 
M, with ^{3^ = g^^^ = M and G^l = 5^2^ = M* ■ The 
parameter 62 retains its bare value since (^2 cos 2ip) ~ 
across the Fermi surface. Introducing the scattering 



amplitudes Ti = 7r7Vi?nimpu,j , 
consistent equations 



?, we obtain a set of self- 



7 i>„ni + ill2 
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where i/„ = w„ + i(5o,7 = Tq -|- F^, M is the (real) SDW 
order parameter affected by impurities, and Asdw is the 
SDW coupling constant. The Neel transition tempera- 
ture Tat to a pure SDW phase at perfect nesting and 



zero disorder is related to Asdw by 1/Asdw = ln(A/T/v_o)j 
where A is the high-energy cutoff. We emphasize that 
both intra-band and inter-band impurity scattering af- 
fect Tn and the order parameter in the SDW channel. 
We also introduce 
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n2 = 



dip cos2(^ 
2^ Vt 



(6) 



with il = W AP + (z>„ + 162 cos 2lp) . Note that 11^ by 
itself depends on the renormalized variables, i.e. Eqs. ^ 
are non-linear self-consistent equations. 

The solution of the set ([5|) gives the Green's function of 
a dirty SDW magnet, which then acts as a bare Green's 
function for the SC system. Since we are only interested 
in Tc, we restrict our analysis to the linearized SC gap 
equation in the presence of impurities and a non-zero 
SDW order parameter: 
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where ( = Fq — F^, 1/Asc = ln(A/Tc^o) is the coupling 
constant in the SC channel, Tc.o is the superconduct- 
ing transition temperature without SDW and at perfect 
nesting, and 
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(8) 
Equation ^ reduces to the gap equation of an ordinary 



dirty 



superconductor if we set M 



0, 



M = 0. Aher- 
we recover the 



natively, at perfect nesting, Sq = I 
resuhs of Ref J^. 

Results and comparison to experiments. We first 

consider the pure SDW state. In Fig. |3] we show that 
both the SDW transition temperature Tn and the or- 
der parameter M are reduced in the presence of impuri- 
ties. This behavior is entirely expected, since both intra- 
band and inter-band impurity scattering are detrimental 
to SDW. A less obvious result is that impurities also af- 
fect the character of the SDW transition at low T. In 
the clean case, the SDW transition is first-order at low 
enough TJ^ Impurities add additional scattering and ef- 
fectively shift T — !■ T + 7, extending the range of the 
second-order transition to smaller temperatures T. Once 
7 exceeds a critical value jcr, the second-order transition 
line extends down to T = 0. The critical jcr is obtained 
in a straightforward way by expanding the last equation 
in ([5]) to order M^ and verifying when the cubic coeffi- 
cient changes sign. For 82 = 0, we obtain analytically 
7cr/ (27rr/v^o) ~ 0.08, where Tn,o is the SDW transition 
temperature at perfect nesting. 

We now use the SDW results as input and solve Eq. ([7]) 
for Tc- To verify whether Tc is reduced or enhanced with 
increasing disorder, it is sufficient to consider small Fq 
and F^ and evaluate AT, = Thirty _2x:ican ^^ gj.g^ ^^.^g^. ^^ 



Ti/Tis!fl. The computations are tedious but straightfor- 
ward, so we skip the details and present our resuhs. The 
phase diagram in the presence of impurity scattering is 
shown in Fig. [T]for on-site impurity potential (Fq = F^) 
and two ratios of T/v,o/T'c,o- We clearly see two differ- 
ent types of behavior in the coexistence phase: Tc either 
increases when impurities are added, or decreases at a 
slow rate. This non-universal behavior can be under- 
stood qualitatively: for a small T/v,o/2^c,o (Fig- Ha)) the 
effects of disorder on T/v and T^ are comparable, and the 
feedback on T^ from the reduction of T/v overshadows 
the direct pair-breaking effect on Tc- For a larger ratio 
Tn,o/Tc,o, the effect of disorder on SDW gets relatively 
weaker, and the direct pair-breaking effect on Tc prevails 
(more specifically, we estimate that Tc increases with dis- 
order when the ratio Tn^q/ (Fq + Ttt) is smaller than the 
ratio Tco/r^). 

In Fig. [Hwe plot ATc = Tf"^y - Tf'''"' as a function 
of the ratio between intra-band and inter-band impurity 
scattering amplitudes. We consider the two sets of pa- 
rameters of Fig. [T] with Sq in the coexistence region and 
compare them with the case when no SDW is present. In 
the latter, ATc < when Fjr is non-zero. We see that 
in the coexistence region ATc is definitely positive when 
Fq/Ftt is large enough, i.e. Tc increases when impurities 
are added into the coexistence state. As expected, this in- 
crease is the largest when r^r vanishes, since in this limit 
impurities are not pair-breaking, but still suppress M. 
When Fq and F^ are comparable, Tc can either increase 
or decrease, depending on parameters, but even when it 
decreases, the rate of the decrease is smaller than that 
for a pure SC state. 

Our results offer an explanation for the non-monotonic 
behavior of ATc = T^"*-^ - Tf"""" as function of doping 
observed in Refs. pj||ll| by adding Zn to LaFeAs(0i_2;Fa;) 
and Cu to Ba(Fei_a;Coa;)2As2, respectively. They also of- 
fer an explanation for the observation in Ref. [I3 that the 
addition of Zn to underdoped Ba(Fei_2;Coa;)2As2 leads to 
a decrease of ATc, but at a slower rate than in the over- 
doped region. For this material, the fact that Zn sub- 
stitution leads to ATc < while Cu substitution leads 



to ATc > may be due to different disorder potentials 
associated with each dopant, leading to different ratios 
Fq/F^. It is also possible that the foreign element not 
only acts as an impurity but also changes the electronic 
chemical potential and/or the Fermi surface geometry. 

We caution that the microscopic coexistence of SC and 
SDW orders has been well-established for 122 compounds 
like Ba(Fei_^Co2;)2As2; but this issue has not been set- 
tled for the 1111 systems like LaFeAs(0i_3:Fa;). Yet, 
even if SDW and SC phase-separate and occupy different 
parts of the sample, we expect some of the physics de- 
scribed here to hold, i.e. that impurity scattering on the 
one hand is pair-breaking and on the other hand tends 
to increase Tc by suppressing the competing SDW phase. 

Conclusions. In summary, we showed that the dif- 
ferent behaviors of Tc with impurity scattering observed 
in overdoped and underdoped iron-based superconduc- 
tors can be understood within the s"' scenario for su- 
perconductivity. While in the overdoped regime Tc is 
quickly reduced with increasing impurity scattering, in 
the underdoped regime there are two competing effects: 
the direct pair-breaking by impurities, which reduces T^ 
and the suppression of the coexisting SDW order param- 
eter, which increases Tc. We demonstrated that, due to 
competition between these two effects, Tc in the coex- 
istence region either drops at a smaller rate or even in- 
creases with increasing impurity concentration, in agree- 
ment with the experimental data. We view this agree- 
ment as an evidence that the gap symmetry in the iron 
pnictides is indeed s'^ . 
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